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^ ' Abstract 

> : 

^ ■ We propose that large quantum fluctuations of the conformal factor drastically 

modify classical general relativity at cosmological distance scales, resulting in a scale 
invariant phase of quantum gravity in the far infrared. We derive scaling relations 
for the partition function and physical observables in this conformal phase, and 
suggest quantitative tests of these relations in numerical simulations of simplicial 
four geometries with 5* 4 topology. In particular, we predict the form of the critical 
curve in the coupling constant plane, and determine the scaling of the Newtonian 
coupling with volume which permits a sensible continuum limit. The existing nu- 
merical results already provide some evidence of this new conformal invariant phase 
of quantum gravity. 



At accessible distance scales, from centimeters to light years, gravitational phe- 
nomena are described quite accurately by the classical Einstein theory. It is very 
well known that this theory is beset with ultraviolet divergences at the quantum 
level and is virtually certain to require drastic modification at the Planck scale. 
A corollary of this severe behavior at ultra-short distances would seem to be mild 
behavior in the infrared. One does not normally think of quantum fluctuations of 
the metric field as important at large distances, and indeed in perturbation theory 
around flat space there is no sign of infrared problems. However, gravitation is a 
long range force characterized by massless excitations that cannot be shielded. If 
the classical spacetime is curved on some characteristic distance scale £, fluctua- 
tions with wavelengths longer than I need not remain small, a fact long ago pointed 
out in Newtonian theory by Jeans. Infrared behavior of fluctuations depends very 
much on the background geometry, and may force large modifications of the classical 
background. 

The quantum analog of the classical background is the ground state. Prior to 
and independent of the question of how to tame the ultraviolet divergences bedev- 
iling the quantum theory, the question of what is its correct ground state presents 
itself. This issue is one that intrinsically involves infrared properties, where we may 
hope to say something sensible in the language of low energy effective lagrangians, 
without knowledge of Planck scale phenomena. 

In order to construct a low energy effective lagrangian for gravity that incorpo- 
rates infrared fluctuations correctly, one must decide first of all what is the relevant 
order parameter field at large distance scales. Here observational cosmology comes 
to our aid to suggest that the FRW scale factor, or more generally, the conformal 
factor of the metric tensor should play an important role. In the classical Einstein 
theory this scalar part of the metric does not propagate. It is determined in terms 
of the matter sources and has no independent dynamics of its own. At the quantum 
level there is a trace anomaly in the energy-momentum tensor of conformally cou- 
pled matter fields. The existence of the conformal anomaly means that the classical 
constraints which fix the scalar part of the metric fluctuations in terms of matter 
sources cannot be maintained upon quantization. In other words, the conformal 
factor becomes unconstrained in the full quantum theory. 

The low energy effective lagrangian must be modified accordingly to take ac- 
count of the trace anomaly and fluctuations of the conformal factor, and this mod- 
ified theory reanalyzed to discover the correct infrared behavior of the quantum 
theory of gravity. In two dimensional quantum gravity, otherwise known as non- 
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critical string theory, the effective action induced by the trace anomaly modifies 
the dynamics of the conformal factor at all distance scales in dramatic fashion. 
Certainly nothing like the fractal behavior and scaling relations of random surfaces 
in 2D gravity can arise from the classical two dimensional Einstein-Hilbert action 
which (being a topological invariant) yields no dynamics whatsoever [1-2] . 

In an earlier paper we obtained the effective Wess-Zumino action induced by 
the trace anomaly of conformal matter in four dimensions. We analyzed this con- 
tinuum effective theory showing that it possesses a non-trivial, infrared stable fixed 
point, characterized by certain anomalous scaling relations [3]. We argued that this 
scale invariant phase is the ground state of 4D quantum gravity approached at dis- 
tance scales much larger than the horizon length of any given classical background. 
Conformal symmetry, apparently broken by the trace anomaly is restored dynami- 
cally by the large fluctuations in the conformal factor at these scales. An important 
consequence of restoration of scale invariance at large distances is the screening of 
the effective cosmological term as measured by the average scalar curvature in the 
ground state. In particular, the two-point correlator of Ricci scalars (R(x)R(x')} 
falls to zero with a certain universal power of the invariant distance between x and 
x' that depends only on the total number of massless fields in the theory through 
the effective central charge of the theory. 

In principle, it should be possible to verify the existence of such a scale in- 
variant phase of 4D quantum gravity by appropriate numerical simulations. In 
this letter we suggest that a comparison of the continuum predictions can be made 
with simplicial simulations of four-geometries with the topology of S" 4 [4-5]. This 
topology includes in particular the physically interesting case of Euclidean de Sitter 
space. We first derive the scaling behavior of the partition function for this topol- 
ogy, or equivalently the behavior of the fixed volume partition function at large 
volumes V. The continuum theory predicts the existence of a critical curve for the 
cosmological term A as a function of the Newtonian coupling k. This implies that 
the corresponding parameters of the dynamical triangulation approach must obey 
a quadratic relation in the scaling limit with the slope of the linear term completely 
determined in terms of pure numbers. 

Moreover in order to obtain the continuum limit corresponding to the scale 
invariant phase of quantum gravity, the Newtonian coupling of the lattice theory 
should be scaled to infinity like \/V for large V. An indication that this is a proper 
infinite volume limit to take is that apparently only in this limit does the average 
scalar curvature vanish on the lattice [4] . A clear and non-trivial test of the existence 
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of the scale invariant phase is that the entropy exponent must be independent of 
the rescaled Newtonian parameter. If this turns out to be the case in the numerical 
simulations, we will have at hand a powerful computational tool for 4D quantum 
gravity in the far infrared, and should be able to measure the graviton contribution 
to the effective central charge of the conformal theory. 

Let us first recapitulate our principal results in the continuum [3,6]. We began 
with the conformal decomposition of the metric g a b(x) = e 2a ~( x )g a i,(x), with g a b( x ) 
a fixed fiducial metric. By consideration of the general form of the trace anomaly 
for conformal fields in four dimensions and taking into account the Wess-Zumino 
integrability condition, we determined the general form of the effective action whose 
a variation is the trace anomaly. Treating this effective action as the fundamental 
quantum action for the a field at large distances and requiring that general covari- 
ance be exactly preserved in the vacuum state of this a theory, we found that the 
total trace anomaly of the full theory must vanish. In other words the absence of 
diffeomorphism anomalies in quantum gravity requires that the vacuum is a scale in- 
variant conformal fixed point where the beta functions of all renormalized couplings 
are zero. 

At the fixed point the effective Euclidean action for a reads: 

Seff = J ^^{^[aA^+I^-ln^^-fe^fa^V^ + fl+Ae 4 -} 

(1) 

where G is the Gauss-Bonnet integrand whose integral is the Euler number, 

X E = ^j d 4 x^ G , (2) 
and A 4 is the Weyl covariant fourth order operator acting upon scalars: 

A 4 = □ 2 + 2i? a6 V a V 6 - fi?D + i(V a i?)V a . (3) 

The quantity Q 2 plays the role of the central charge at the infrared fixed point and 
is proportional to the coefficient of the Gauss-Bonnet term in the quantum trace 
anomaly. 

The physical metric at the conformal fixed point becomes 

g ab (x) = e 2a ^g ab (x), (4) 

with a determined by the condition that the Einstein-Hilbert action have its canon- 
ical scale dimension in this metric. This condition gives a quadratic equation which 
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fixes a in terms of the central charge: 



a = 



1-Jl- 



2 



(5) 



With a determined in this way the condition that the volume (cosmological) term 
have dimension 4 requires that the cosmological and Newtonian couplings satisfy 
the relation 1 , 



A 



K 



f(Q 2 ) 



a 



HQ') = q2 



4a 2 6a 4 
+ ~Q 2 + ~W 



(6) 



All of these results were obtained in the continuum by treating the metric g a b 
as fixed. In other words the transverse, traceless sector of the theory containing the 
physical spin-2 gravitons was neglected completely. Our basic hypothesis is that 
these relations remain true in the infrared when the graviton modes are included, 
up to a possible renormalization of the value of Q 2 . More precisely we assume that 
integration over the transverse graviton modes generates an effective action for a 
which, when expanded in powers of derivatives, has the same form as (1) but with 
renormalized coefficients. This assumption we call "infrared conformal dominance." 

This is not at all unreasonable from a Wilsonian effective action point of view. 
Consider the functional integration over transverse gravitons (in other words over 
conformal equivalence classes of metrics), as well as matter fields, with both in- 
frared and ultraviolet cut-offs, i and a. At short distances, graviton effects may 
grow uncontrollably due to the presence of the dimensionful Newtonian coupling 
k, so that a cannot be taken to zero in the effective action. Conversely, at large 
distance scales, the transverse, tracefree fluctuations should be expected to become 
less important, so that the effective action should be regular as the infrared cutoff 
I is removed. If this is the case, then an infrared stable renormalization group 

fixed point of the effective low energy theory is approached as t — > oo. Scale 
invariance at this fixed point then requires that the low energy effective action must 
be of the form (1) when expanded up to four derivatives of a. The uniqueness of 
the effective action (1) at the infrared fixed point is still not sufficient to guarantee 
that the relations (5) and (6) hold after the integration over the graviton modes. 
For these relations to be unmodified in terms of Q 2 one needs to assume that the 



1 Note a typographical error in eq. (3.17) of ref. [3] and the difference in the 
definition of A. 
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infrared fluctuations of the gravitons are subdominant compared to the fluctuations 
in the confer mal factor. 

In (1), we have dropped a possible C 2 hcd a term. This is justified by scale 
invariance at the fixed point, which requires the vanishing of the beta function of 
the Weyl-squared coupling, so that the coefficient of the C 2 bcd term in the trace 
anomaly must vanish identically at the fixed point. The same reasoning eliminates 
the coefficient of a local R 2 term in (1), which can be checked explicitly at the level 
of the a theory alone. Then the only unknown parameter in the continuum theory 
is the contribution of gravitons to Q 2 . We have calculated this contribution to Q 2 
in perturbation theory in both the Einstein and Weyl theories and found values 
around 8, which lead to an a ~ 1.2 [6]. 

In order to derive the scaling behavior of the partition function of the effective 
a theory, subject the a field to the constant shift [2] 

a^a+- (7) 

a 

and use the translational invariance of the integration measure [Da] to find: 

Z(k, A) = J [Va]e- S ^f [(7] = e-^- x ^Z(Ke- 2uj ,Xe 4uj ) . (8) 

In what follows we restrict to the topology S 4 for which x B = 2. It is convenient 
to define also the partition function at fixed volume 



Z(k; V) = J [Do]e xv - s °"W S ( J d 4 x^e 4aa - V^j 
Then performing the translation (7) above, we obtain 2 

Z(k; V) = e- 2 ^ +2 ^Z(Ke- 2 "; e~^V) 



(9) 



(10) 



where in the second line we put e 4ai oc V. 

From these continuum results we turn now to the numerical simulations. The 
numerical method that has proven most fruitful up until now is that of "dynamical 



2 Scaling relations of a similar kind were derived in ref. [7] for conformally self- 
dual metrics. However, the physical meaning of the second operator identified in 
this work with the "volume" is unclear to us. The introduction of this new operator 
is also the reason why the critical relation (6) was not obtained. 
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triangulation," a variation of Regge calculus in which geometries are constructed by 
gluing together fundamental simplices of fixed volume [8] . The four-simplices share 
common faces with their neighbors which are D — 1 = 3 dimensional simplices, i.e. 
regular tetrahedra of edge length a. The angle between two tetrahedra faces sharing 
a triangle is 

9 = arccos(-^) = 1.3181161 (11) 
in D = 4 dimensions. The volume of a fundamental J-simplex is 

vj=aJaj= j<\fW' <i2> 

where a is the lattice spacing, so that the total volume of the simplicial manifold is 

Jd 4 x^^N 4 V 4 (13) 

if -/V4 is the total number of 4-simplices in the configuration. 

As in the Regge approach, space is regarded as flat inside the D = 4 simplices 
with all the curvature residing on the D — 2 = 2 dimensional hinges, i.e. equilateral 
triangles with the same fixed edge length. If rii is the number of 4-simplices sharing 
a given equilateral triangle i, then the deficit angle Si is given by 

Si = 2tt - mO , (14) 

and the Einstein-Hilbert action takes the value, 

/ dtx y/tj R^J2 SiV2 = ( 2?riV2 " 10dN *) V ^ ( 15 ) 
"* i 

where J2i n i = 10A^ 4 has been used (since each 4-simplex has 10 triangles in its 
boundary) . 

Dynamics is now specified by giving an action for each simplicial triangulation 
T of the form S(T) = J2j kjNj(T), with each triangulation otherwise having equal 
statistical weight. Actually not all of the Nj are independent. In order for the sim- 
plicial complex to approximate a continuous manifold the numbers of J-subsimplices 
(for J = 0, 1, . . . , D) must satisfy some relations, called Dehn-Sommerville relations 
[9]. In addition we have the Euler relation, 

Y j {-) J N j = Xe . (16) 

j=o 
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The net result is that only two of the Nj are independent and the action may be 
taken to be of the form, 



S(T) 



k 2 N 2 (T) + k 4 N 4 (T) . 



(17) 



Comparing this with the Einstein-Hilbert action, and using the substitutions (13) 
and (15) leads to the following identification of the simplicial action parameters 
with those of the continuum, 



where the numerical coefficients come from eq. (12). 

Since we are interested in the continuum infrared fixed point of the lattice 
theory, we do not add any higher derivative couplings to the action. These should 
correspond to irrelevant operators in the infrared, and in any case the coefficients 
of possible R 2 and Weyl-squared terms in the action vanish at the conformal fixed 
point. The trace anomaly induced action (1) is not to be added to the lattice 
action (17) either. It is nonlocal in the full metric (4) and should be generated 
dynamically by the quantum fluctuations of the simplicial geometries, in analogy 
with the situation in the two dimensional case. 

Because the number of triangulations with fixed S 4 topology which can be 
made from a given number of 4-simplices is exponentially bounded with respect 
to N4 [4-5] , the partition function of the dynamical triangulation approach, 



Z DT (k 2 , h) = e" S(r) = E Z ^ N ^~ k4N4 & E e-^-WW* (19) 



must exist in a region of the coupling constant plane k± > k%{k 2 ). By approaching 
the boundary of this region from above one can hope to arrive at a continuum 
limit in which physical correlation lengths go to infinity when expressed in lattice 
units. In other words, one is searching for a critical curve in the {k 2l k^) plane 
corresponding to a second order phase transition of the lattice theory defined by 
(17) and (19). 

The first implication of the continuum a theory for the lattice simulations is 
that at the infrared fixed point one has the relation (6) 



v^tt a 2 



4 K 

a 2 
4 K 




(18) 



r 
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which determines the critical curve in the (k 2 , k 4 ) plane, 

kt(k 2 ) = —k 2 + Vbf(Q 2 )k 2 2 , (20) 

7T 

where use of eqs. (18) has been made. 

This relation should hold for the lattice parameters in the infinite volume (con- 
tinuum) limit. In any finite volume simulation there will be an additive renor- 
malization of the cosmological term which scales to zero with the lattice length a. 
Hence the intercept of the curve k 4 (k 2 ) will not vanish in general, for finite volume. 
However, the slope of the linear term of the critical curve (20) is = 2.0978469, a 
pure number independent of Q 2 . Actual simulations with N 4 ~ 10 4 seem to indi- 
cate a critical curve k 4 (k 2 ) which is approximately linear with a slope slightly more 
than 2 [10]. The relation (20) has a quadratic term as well which could be used to 
determine Q 2 , in principle if the simulations are run with high statistics and large 
volumes. However since Q 2 is unknown, and possibly large, and f(Q 2 ) — > for 
large Q 2 , the quadratic term could be difficult to measure. 

A quite different prediction of infrared conformal dominance and better way 
to measure Q 2 is provided by the finite volume scaling relation (10). Translating 
this continuum relation to the lattice by using (18) we obtain the following scaling 
relation for the fixed volume partition function at large volumes, 

Z(k 2 - N 4 ) ~ Z(k 2 )Nr 3 e k ^ N * , (21) 

where 

k 2 = k 2 VN 4 , (22) 

and 

7(Q 2 ) = 2 - ^ • (23) 

The content of the scaling relation (21) is that, when k 2 is scaled to zero with the 
square root of the volume as in (22) keeping k 2 fixed, the prefactor Z becomes 
volume independent and the entropy exponent 7 depends only on the effective 
central charge Q. Therefore, a clear test of this scaling relation is that 7 must be 
independent of the parameter k 2 , if our hypothesis of infrared conformal dominance 
is correct. Otherwise, 7 would acquire non-trivial k 2 dependence. If 7 is indeed 
independent of k 2 , a measurement of 7 would provide a non-perturbative way to 
compute the graviton contribution to the central charge. Note also that the entropy 
exponent goes to —00 in the semiclassical limit Q 2 — > 00, while 7=1 for Q 2 = 
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4. In the latter case one would expect logarithmic behavior in analogy with the 
c = 1 case in two dimensions. For Q 2 < 4 the exponent a in (5) becomes complex 
and the theory could exhibit a phase transition with qualitatively new phenomena. 
Perturbative calculations of the graviton contributions [6] lead to the value 7 ~ 
-1.3 . 

The simplest observable with nice scaling behavior is the average curvature (R) 
defined by [4]: 



/ d A x^fgR\ _ 2ttV 2 



' ° lnZ(k 2 ;N 4 )-- 



N,\ dk 2 " 7r 



(24) 



where the last equality holds for fixed volume and we used the relations (13), (15) 
and (17). Now inserting the scaling behavior (21) and the expression (20) one 
obtains: 



1 



2v / 5/(Q 2 )fc 2 + -l r 4UnZ(fc 2 ) 
V N 4 dk 2 







(25) 



where the proportionality factor depends only on Q 2 and the rescaled k 2 . This shows 
that the average curvature scales to zero with the square root of the volume. Hence 
scaling k 2 — > with k 2 fixed yields a large volume continuum limit consistent with 
naive dimensional analysis. This removes the main obstruction to interpretation 
of the numerical simulations for the continuum theory mentioned in ref. [4], since 
(i?) lattice = a2 (-R)continuum vanishes as k 2 ~ N A — > 0. This is consistent with the 
linear behavior of (R) with k 2 observed numerically in ref. [4]. 

Another observable that has been used to search for a continuum limit at a 
critical value of k 2 is the integrated curvature-curvature correlator, 



1 & 



2 



Vdk 2 



In Z(k 2 ;N 4 ) = J d 4 x (R(x)R(0)) - V(R) 2 . (26) 



If we substitute (21) into this expression and use (20) and (22), thereby scaling 

_i 

k 2 — > like N 4 2 , this quantity goes to a constant in the large volume limit, in 
contrast to the divergent behavior at a finite value of k 2 that has been suggested in 
the searches for an ultraviolet fixed point. 

Finally we would like to mention that considering only the first term of the 
curvature-curvature correlator in the r.h.s. of (26) is also interesting from our 
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present perspective. This is 




iV 4 ) Z(k 2 ;N 4 ) . 



(27) 



If we substitute (20), (21) and (22) into this last expression, a short calculation 
shows this correlation function also goes to a finite constant, independent of N4 
in the large volume limit. This means that the correlator (R(x)R(0)) must fall off 
faster than \x\ for large \x\, in order for the integral over x to converge. Contrast 
this convergent behavior of (27) in the conformal invariant phase of 4D gravity to 
the behavior of the same quantity in the classical theory. In that case the trace of 
the classical Einstein equations, R = 4A rigidly fixes the scalar curvature in terms 
of the cosmological constant. Hence (27) must diverge linearly with the volume N4 
in the classical theory, if the cosmological term is non-vanishing. 

We conclude that if the quantum fluctuations of the conformal factor described 
in this letter really dominate in the far infrared, then any cosmological term in the 
continuum theory is screened completely at the largest distance scales by these fluc- 
tuations, or in other words, the effective cosmological constant in the scale invariant 
vacuum of 4D quantum gravity must vanish identically. 
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